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We show that given three hermitian matrices, what one could call a fuzzy repre- 
sentation of a membrane, there is a well defined procedure to define a set of oriented 
Riemann surfaces embedded in M 3 using an index function defined for points in M 3 
that is constructed from the three matrices and the point. The set of surfaces is co- 
variant under rotations, dilatations and translation operations on M 3 , it is additive 
on direct sums and the orientation of the surfaces is reversed by complex conjugation 
of the matrices. The index we build is closely related to the Hanany-Witten effect. 
We also show that the surfaces carry information of a line bundle with connection 
on them. We discuss applications of these ideas to the study of holographic matrix 
models and black hole dynamics. 
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I. INTRODUCTION 

The discovery of D-branes [Tj introduced a huge class of new geometric objects in string 
theory. It was quickly realized that the coordinate positions of these geometric objects are 
matrices, and hence that their positions can become smeared in the same way that typical 
wavefunctions in quantum mechanics do not have a well defined position and momentum 
due to the uncertainty principle. This fuzzyness of the D-brane position occurs when the 
matrices that describe the positions of D-branes do not commute. When the said matrices 
commute, the positions of the D-branes can be identified with the eigenvalues of the matrices 
themselves. 

A second route to obtaining noncommutative coordinates for branes arises from the light- 
cone quantization of the membrane [2J. Indeed, just the lightcone description in the classical 
theory itself leads to such a prescription. In that case the coordinates that don't commute 
describe the internal coordinates of the membrane itself. The idea of Goldstone and Hoppe 
is that in the lightcone quantization, the supermembrane acquires a non-degenerate Poisson 
bracket on its spatial worldvolume. These Poisson brackets on the membrane coordinates 
are then approximated by commutators of finite matrices, so the matrices become the inter- 
nal coordinates of the membrane itself. This is a UV truncation in the degrees of freedom, 
so it regularizes the membrane theory on the lightcone. Indeed, the supermembrane version 
of this construction j3] was one of the pieces of evidence that was given in the construction 
of the BFSS matrix model describing the full M-theory on the lightcone quantization [3]. 
The geometric object we call the membrane is supposed to appear in the limit where the 
size of the matrices goes to infinity, with small commutators, but it is unclear if a geometric 
object that is a membrane can be defined for finite matrices or not, and whether it is a sharp 
geometric object or a very fuzzy object. 

All of the considerations above occur at the classical level. There is no need to invoke 
quantum mechanics to have these noncommutative geometric effects occur. We can then 
ask the following question: given a collection of matrices that don't commute, is it possi- 
ble to construct a set of surfaces associated to them that would represent the membrane 
worldvolume geometry (or multiple membranes) embedded in flat space without taking an 
infinite size matrix limit? 

We will answer this question in the affirmative for the case where we are given three 
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hermitian matrices X, Y, Z that generically do not commute. These would then represent 
an embedding of a membrane in M. 3 . We will also show that this easily generalizes to 
embeddings on a plane wave in the lightcone where the plane wave has a transverse M. 3 set 
of coordinates. Our construction produces a collection of closed oriented surfaces embedded 
in IR 3 . Indeed, we will show that these membranes not only carry an orientation, but that we 
can also deduce that they carry vector bundles on them and hence behave like D2-branes. We 
will also show that the topology of the membranes is continuous when X, Y, Z are deformed, 
so topology transitions are not instantaneous and require going through singular geometries. 
This implies that the associated brane charge is conserved. 

The main motivation to understand this problem in detail arises from the study of sim- 
ulations of black hole formation in matrix models performed in [5]. The data obtained 
there at the end of the evolution is exactly of the sort above: a collection of (somewhat 
random) matrices that do not commute. Although one can try to find the brane positions 
by diagonalizing each matrix, the end result is not amenable to easy visualization in higher 
dimensions. After all, what are we supposed to do with the noncommutative information? 
The approach in [6] is to find some approximate locations for DO branes that minimizes the 
non-diagonal matrix elements given this choice of basis and then throws the off-diagonal 
information away. Then it uses those DO branes as a proxy for the geometric object. This is 
a very nice idea. However, this type of description does not describe the orientation of the 
extended membranes nor can it be used to determine the topology of the brane configuration 
except in the large matrix limit. This is especially hard if one is near a topology change. It 
also seems to indicate that the result is very fuzzy and the topology of the branes is in the 
end given by the topology of a set of points. 

Our approach to this problem is very different. We start from the BFSS jl] and BMN 
matrix models [7] and to simplify matters in the discussion, we orbifold the problem suffi- 
ciently so that in the end we can deal with a reduced model where only three matrices are 
required, so instead of starting with a system with 16 supersymmetries, we go to a system 
with 4 supersymmetries. We do this by taking a model which is the dimensional reduction 
of Zfc supersymmetric orbifolds in four dimensions that are chiral. The model we need is 
then obtained by taking a quiver where only fractional branes on one of the nodes of the 
quiver are present, this is, we concentrate on the reduction of pure N — 1 SYM reduced to 
matrices (there are 3 of them that are dynamical, which we call X, Y, Z). To explore the 
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geometry we then add a fractional brane probe in one of the nodes of the quiver diagram that 
intersect the node where our matrices are located. The main idea is to ask what the D-brane 
probe sees given a generic set of coordinates as given by the X, Y, Z matrices. Within the 
dynamics of the matrix model orbifold, there are either bosons connecting the probe brane 
to the configuration, or fermions. The interesting degrees of freedom to define the geometry 
end up being the fermions. Indeed, one can show that the spectrum of fermions connecting 
the probe to the fuzzy object can be obtained by diagonalizing a simple hermitian matrix 
obtained from the X, Y, Z and the coordinates of the probe. The matrix can be thought of 
as an effective Hamiltonian and it is given by 

H eff = (X-x)a x + (Y-y)<j y + (Z-z)a z (1) 

A technical point is that the eigenvalues of H e jf are n °t all positive nor all negative. In 
second quantizing the fermions the positive eigenvalues are associated to raising operators, 
and the negative eigenvalues to lowering operators. The conjugate modes come from the 
anti-chiral fermions. The absolute value of the eigenvalues then serves as the mass for the 
fermions. This mass can be thought of as the length of a fermionic string connecting the 
probe to the configuration, so it gives a notion of distance, and the minimal eigenvalue is 
the shortest distance to the configuration. 

The important technical point of this paper is that the eigenvalues of H e ff can cross 
zero and this depends on the position of the probe. The geometric surface locus of the 
matrix configuration is described exactly by the locations where one of the eigenvalues of 
H e ff vanishes. Counting the number of positive eigenvalues of H e ff versus the negative 
eigenvalues can indicate the number of such crossings of zero, and this can be used to define 
an index I(x, y, z)x,y,z- The plane M. 3 represented by the coordinates x, y, z is then colored 
by the index of the location. The index is locally constant and can only change if one of 
the eigenvalues crosses zero. Any path connecting two points with different index must have 
zero crossings: it is impossible to avoid them by taking a clever path. This indicates that 
the surfaces obtained this way are closed. At any crossing, we also get an orientation: from 
higher index to lower index, so the surfaces are oriented. 

A second point that is worth mentioning is that at the zero crossing a raising operator 
becomes a lowering operator, or viceversa. Thus if one follows a vacuum of the fermion de- 
grees of freedom on one side of the brane continuously to the other side we get an anomalous 
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creation of fermionic strings. This is a generalization of the Hanany-Witten effect [8] and 
it represents anomalous creation of branes by branes (for more details on the relations to 
topology and anomaly inflow of this effect see [5]). 



The paper is organized as follows. In section |TT] we review the BFSS and BMN matrix 
models and some of their orbifolds. This is used to find an effective Hamiltonian for chiral 
fermions in the presence of a configuration of three matrices plus a probe, which was already 
written in equation Q. Next, in section III we show how the number of positive versus 
negative eigenvalues of the effective hamiltonian can be used to define an index function 
given a position of a probe. The index vanishes when the probe is at infinity. The locus 
where the index changes defines a collection of surfaces. We show various properties of the 



index. In section IV we show that fuzzy spheres give configurations where on some loci the 
index is non-zero. The associated surfaces are spheres. We also show how to construct torus 
embeddings in IR 3 . Finally, we show that the surfaces carry the information of a line bundle 
on them. This makes the behave like D-branes. In section [V] we show how to generalize the 
index to a linking number between two such configurations. The linking number ends up 
counting the number of strings that are created by trying to separate the two configurations 
in a generalization of the Hanany-Witten effect. Equally, the index counts the number of 
strings that are created when bringing the configurations together from infinity. We show 
that in a special case of fuzzy spheres for the BMN model, that the data on when the fermion 
zero modes appear provides additional evidence that the surfaces carry a line bundle on them 
and behave like D2-branes. We apply these ideas to study matrices obtained by numerical 



studies of the BMN matrix model in section VI We give applications to understand the 
polarization of black holes into membranes and we also give applications where the Hanany- 
Witten effect can stop probe D-branes: we show that as N increases, the Hanany-Witten 
effect strings end up storing more energy parametrically in N than the probe, so they are 
enough to show that the black hole will stop the probe inside it and it will not come out at 



the other side. Finally, in section VII we conclude. We also have appendices where some of 
our more elaborate computations and conventions are described. 
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II. THE BFSS AND BMN MATRIX MODELS, AND THEIR ORBIFOLDS 

The BFSS matrix model [I] is simply the dimensional reduction of 9 + 1 SYM to + 1 
dimensions. It is a gauged matrix quantum mechanics of 9 adjoint matrices <p l with a gauge 
symmetry by similarity transformations. An amazing aspect of this dynamical theory is 
that it is a description of M-theory in the discrete lightcone quantization and in the limit 
iV — > oo and in an appropriate double scaling limit it gives a quantization of M-theory in 
flat space in a lightcone quantization. The action is given by 



Sbfss = I dt Tr 



9 1 

+ ^-(2i?)(*V[^>]) 



j,k=i 

(2) 



3=1 

As written, the action depends on a parameter R. This parameter can be eliminated by 
a rescaling of the variables and time and it can be replaced by h. Thus the BFSS matrix 
model itself has no intrinsic scale at the classical level. Indeed, one can check that even h 
can be removed because the action has a classical scaling symmetry. 

The simplest version of the model is for 1 x 1 matrices. In this case, the classical config- 
urations are given by a point in IR 9 with a velocity, and the fermions just add degeneracy to 
these states. This degeneracy gives the correct degrees of freedom for a graviton supermul- 
tiplet in 11D [I]. The M. 9 describes the transverse directions to the lightcone, and the rank 
of the matrices is the amount of lightcone momentum. Such a configuration is a DO-brane. 
Ground states in the classical theory in general correspond to configurations of commuting 
matrices, where N such DO-branes are located on M. 9 . 

One of the important things about the BFSS matrix model is that it is capable of describ- 
ing extended objects. Indeed, one can describe D2-branes and higher order D-branes from 
these configurations. These are easy to see in the infinite N limit [10J, as central charges 
can be activated in the supersymmetry algebra that encode such objects of infinite extent. 
Such solutions lead to effective noncommutative field theories. A modern introduction to 
the geometric interpretation of these developments can be found in [llj. 

One can also check that matrix configurations source the various supergravity fields at 
long distances, and that the couplings to weakly curved backgrounds give us a way to com- 
pute the currents and the multipoles with respect to the brane charges of the configurations. 
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This was very systematically developed in the works of Taylor and collaborators [12J. A 
review of the BFSS matrix model where all of this is very clearly addressed is in |13j . 

Finite matrix configurations can also behave like extended D-branes. The simplest ex- 
ample of such configurations are fuzzy spheres [14] , where three of the matrices are made 
proportional to angular momentum matrices. These have been studied extensively. An im- 
portant question to ask is if these geometries survive at finite N, or if they are only well 
defined strictly when N — > oo. We will show in this paper that precise geometries can 
be described even for finite N. However, since in the BFSS matrix model in principle one 
can also describe all other D-branes in type IIA string theory, a random configuration of 
matrices would be too complicated: it would probably encode somewhat random extended 
D-branes of type IIA theory. It would be nice if we could reduce the problem to just studying 
surfaces in three dimensions, where only three of the <f> matrices matter, and the other six 
are eliminated somehow. 

A very simple way of doing this is by realizing that the BFSS matrix model can also be 
thought of as the dimensional reduction of N = 4 SYM in four dimensions down to + 1 
dimensions. If we manage to reduce the supersymmetry from N = A SYM to just N — 1 
SYM, then instead of having 9 matrices l , we would get only three matrices, those that 
arise from the dimensional reduction of the gauge field connection. In that situation, the 
DO branes would be confined to an IR 3 , rather than an M 9 , and we might expect that we 
can only describe D2 branes, as any higher dimensional even brane would have too high a 
dimension to fit in three dimensions. 

A simple way to achieve this truncation and to keep a full geometric interpretation of 
the system in terms of string theory is to take a supersymmetric orbifold C 3 /Z^. These 
are described by quiver theories which can be constructed by the techniques developed by 
Douglas and Moore [To]. What matters for us is that we can end up with such theories 
where only the dynamics of N = 1 SYM matter. 

This would be the theory of N identical fractional branes at the orbifold singularity. A 
simple explanation of how those field theories can be built and studied is found in [16] . 
The geometric interpretation in terms of fractional branes and intersection theory of those 
objects can be found in [T7] . 

For simplicity we can choose a action that gives rise to chiral theories where between 
any two nodes in the quiver there is at most one chiral field connecting them: this way if 
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we add a probe for a different fractional brane, we can get a single chiral multiplet worth of 
fields connecting the probe to the configuration. The particular example of an orbifold we 
can choose is given by acting on C 3 given by coordinates «i, «2, a 3 and acting with the 
defined by the identifications ot\ — > uai, &2 — > w 2 «2, «3 — > u~ 3 a3, and u = exp(2ni/k) is 
a primitive root of unity. Many other orbifolds will have similar properties and the precise 
details of the orbifold are not important at this stage. All that we need can be visualized 
by a simple sub quiver diagram 



where the arrow indicates a single chiral multiplet. 

The advantage of having a single chiral field is that the fermions are represented by a two 
component Weyl spinor, and the 7 matrices appearing in the BFSS matrix model reduce to 
the four dimensional gamma matrices for such spinors: those are just the Pauli matrices. 
The details of the reductions are shown in the appendix. The other advantage of having 
chiral fields is that they carry anomalies in four dimensions, thus it is natural to assume that 
they might encode a lot of topological information even in the reduction to + 1 dimensions. 

Upon such a reduction, we end up with fermion terms where we only involve four dimen- 
sional 7-matrices. Indeed, if we reduce to a single chiral multiplet, then we can think of the 
7-matrices themselves as Pauli matrices. The effective action is then 



where if ip is chiral, then ip* is antichiral. Again, R is meaningless as it can be redefined 
away, and the classical symmetries of S or b have the same properties as those for Sbfss- 
The action above is a shorthand: it is the same action of the BFSS matrix model, but the 
matrices are restricted by the orbifold conditions [Po] Q 

1 In practice this can be done keeping the form of the action fixed and adding information about the matrix 
restrictions by using a crossed product algebra |18| . This will produce a set of orthorgonal projectors 
for each node of the quiver, and the commutation relations with these projectors will recover all the 
information of the quiver diagram. For example the traces of the projectors will recover the rank of the 
gauge groups on each node . 




Probe 



(3) 



orb 




(4) 
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The new advantage is that now we only have to deal with three hermitian matrices 1 ' 2 ' 3 
rather than 9. Also, the Pauli matrices are easier to handle than the 9 dimensional gamma 
matrices. 

The question is then if given 1 ' 2,3 hermitian matrices, can we associate a collection of 
D2-branes in a specific geometric configuration in IR 3 to it? To the extent that we can, we 
can then uplift any such intuition to 9 dimensions and understand better how membrane 
geometries arise in the BFSS matrix model. 

Another useful matrix model to consider is the BMN matrix model [7j. That model 
describes M theory on a plane wave in the discrete lightcone quantization. Its action is 
given by 



9 

u mass 




(5) 




\ (9 * f 7l23* 



(6) 



which is a mass deformation of the BFSS matrix model. Again, if we look at 1 x 1 matrices, 
the configuration space is IR 9 , but there are no flat directions: there is a quadratic potential 
in the Hamiltonian. This is as it is supposed to be: it just reflects the fact that there is 
a gravitational potential in the plane wave geometry. This model also can be obtained by 
a dimensional reduction of N = 4 SYM to + 1 dimensions. We need an SU (2) invariant 
reduction on a sphere [TH]. Again, dealing with full 9-dimensional matrices is not very 
intuitive, so we can play the same orbifold trick on </> 4 " 9 to get rid of those matrices. Again, 
we can get rid of R and we can choose units so that ji — 3, but then we are not free to rescale 
h to be whatever we want any longer. Thus the BMN matrix model does have a parameter 
h, also when we orbifold. One can go towards the classical regime, and again, if we give 
three matrices 1 ' 2,3 we can ask: is there a way to associate a geometric D2-brane to such a 
configuration? Indeed, once we take the BMN matrix model at finite N, the ground states 
are made of collections of concentric fuzzy spheres. These are such that the <fi themselves 
are angular momentum matrices with canonical normalization. 

The answer both here and in the BFSS matrix model will be yes: one gets an associated 
set of surfaces in both cases. The surfaces will be slightly different in the BFSS versus the 
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BMN matrix model. The new ingredient in the BMN matrix model is that the fermions 
get a contribution proportional to 7 1 ' 2 ' 3 oc aW 3 oc 1 in their mass. This reflects the fact 
that in the 11 dimensional maximally supersymmetric plane wave background there is a 
non trivial background flux. Such a contribution changes slightly the shape of the branes 
that one associates to the configurations and this is essentially due to the Myers effect [20J : 
branes are polarized in the presence of RR backgrounds. In particular a DO brane becomes 
polarized into a sphere. 

The essence of this article is to look in detail at the fermion degrees of freedom to 
understand the geometry of branes. Our technique is that given three X matrices (these are 
identified with 1 ' 2,3 for one of these models and will be called X collectively, or X, Y, Z if we 
want to name them individually), we will ask: what would a probe point like DO-brane see? 
We ask the question from the point of view of the fermion degrees of freedom that connect 
it to the configuration we are studying. This is encoded in the dynamics of the fermions 
themsleves. The problem reduces to studying an effective Hamiltonian given by 



The last term is for the BMN model and it describes the additional contribution to the 
effective Hamiltonian from flux. The value comes from a choice of orientation of the flux and 
how it relates to the chirality of the fermions (the sign choices when we take 7 1 ' 2 ' 3 ). Solving 
for H e ff gives us the energies of the fermions that connect the object to the configuration. In 
the full dynamics, the wave function solutions for H e ff are second-quantized, as the fermionic 
objects ip become operators when we turn on quantum mechanics. This is important for the 
physical interpretation of the membranes. 



As we have described previously, the geometry in the BFSS and BMN matrix models 
is encoded by matrices of dimension one. In this section we will work exclusively with the 
BFSS matrix model type of dynamics. We will only add some passing remarks at the end. 
To understand how a generic object of the model looks geometrically (a general matrix 
configuration), we can ask the question by adding a point-like probe. This is, we want to 
extend the size of the matrices by one, by taking a direct sum with a zero-brane probe. This 
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III. THE INDEX: ADDING A DO-BRANE PROBE 
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just means that we put the NxN matrix configuration and embed it into an (JV+1) x (JV+1) 
matrix in the upper left corner, we add the eigenvalues in the rightmost bottom corner and 
add zeros everywhere else. This is, we have a new auxiliary configuration where 



In our problem, the matrices X, Y, Z share their properties with the BFSS matrix model: 
they are three hermitian matrices. 

The essence of the geometric characterization of the general matrix will then be encoded 
in the observations of a spectator brane. The spectator brane will only be allowed to ask 
questions related to the dynamics of the matrix models themselves and in particular of the 
degrees of freedom that connect the extra eigenvalue to the matrix configuration. 

In general there are two classes of modes that connect the extra eigenvalue to the con- 
figuration: bosonic degrees of freedom and fermionic degrees of freedom. We will restrict 
ourselves to the fermionic degrees of freedom. The questions we will ask depends on the 
position of the extra eigenvalue probe. 

When we look at the fermions, we decompose them as follows 



where our goal now is to ask what are the energies associated to the modes ip. Notice that 
this picked a very particular component of the fermions and not the other. This can be 
justified completely in orbifold models, as we discussed previously, but orbifolds are not 
really required to make this argument. All we need is the subquiver diagram that enforces 
the restrictions of the matrices defined by equations Q8J) and (J9|. To do this carefully, we are 
choosing the probe to be a different fractional brane than the matrices X represent. The 
chirality of these modes indicates that if the branes were four dimensional fractional branes, 
then they would intersect for sure if we think of fractional branes as higher dimensional 
branes wrapped on collapsed cycles. The intersection properties of the fractional branes 
represent the intersection properties of the collapsed cycles [TF] . Notice also that we did not 
put fermions in the bottom leftmost corner: this is our chirality assumption for the arrow. 

The obvious question to ask first, is if there is a definition of distance from the eigenvalue 
probe to the matrix configuration. The way to ask that question is to look at the spectrum of 






(9) 
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fermions connecting the probe eigenvalue to the matrix configuration. The eigenvalue probe 
is located at x, y, z and the three matrices X, Y, Z are three N x N hermitian matrices. 
Given our three matrices X, Y, Z, this is described by the following effective Hamiltonian: 

H eff = (X - xl N )a x + (Y - yl N )o y + (Z - zl N )a z (10) 

where we are being pedantic in stating that x is multiplying the identity matrix of N x N 
matrices. We will omit this in the future. The structure of how the Pauli matrices appear 
for chiral multiplets is derived in the appendix |Bj 

This is the Hamiltonian of the fermionic degrees of freedom connecting the probe brane 
to the rest of the configuration. The origin of this Hamiltonian is seen from the term in the 
full Hamiltonian given by 

TrWYlXiM (11) 

when evaluated in the configuration X, Y, Z. The Hamiltonian above describes the mass 
term for the off diagonal modes of the fermion if) that are charged under the gauge group 
of the extra eigenvalue probe. The gamma matrices in three dimensions are given by the 
Pauli matrices, whereas the dependence on x and X, etc, comes from direct evaluation of 
the commutators. 

We can think of this as a Hamiltonian in a tensor product space Hilbug = HUbx <8> 
Hilb^i of an N— dimensional Hilbert space times a spin one half object ( a single q-bit). 
H e ff is covariant under unitary transformations of HUbx- This is, we have that under 
U G Aut(Hilb^), we can consider this inducing an automorphism of Hilb^ ig by U ® 1. 
The automorphism takes X UXU' 1 , Y -»■ UYU' 1 , Z UZU' 1 and H eff -»• (U ® 
l)iJ e //(?7 _1 (g> 1) which shows that the spectrum of H e ff is invariant under such rotations. 
This is inherited from the gauge transformations of the original matrix model. What is 
important is that the spectrum of H e ff is gauge invariant. 

As is usual in string theory, the off diagonal modes connecting a sub configuration to 
another are considered to be strings, once they are quantized. The typical energy of a string 
of length I is given by a'£ where a' denotes the string tension. Hence, in our effective 
Hamiltonian, we can denote the distance from the probe brane located at (x, y, z) to the 
configuration by the eigenvalues of the effective Hamiltonian H e ff. The reason to look at 
fermions is that fermionic Hamiltonians do not have tachyons. Thus technically all energies 
are positive and thus the notion of distance is positive. This is also true for string states: 
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open string fermions in the NSR superstring appear in the Ramond sector for open strings. 
The zero point energy of the fields cancels between bosons and fermions on the worldsheet 
(they have the same boundary conditions) and the only contribution to the energy of the 
string is from the classical stretching between the ends of the strings. 

The eigenvalues of H e ff themselves can be positive or negative, so we interpret the pos- 
itive eigenvalues as the frequencies of creation operators, and the negative eigenvalues as 
frequencies of lowering operators once we second-quantize. The absolute value of the spec- 
trum of H e ff is then the list of distances from the probe brane to the object when interpreted 
as strings. Obviously, if we have more than one distance, the object with respect to which we 
are measuring distances should be considered to be an extended object. The minimal eigen- 
value of the spectrum thus obtained should give us the minimal distance to the extended 
configuration. 

Notice that the Hamiltonian i? e // is covariant under rotations and translations. This 
is inherited from the symmetries of the original BFSS Lagrangian. More importantly, the 
Hamiltonian is also covariant under rescalings (if we rescale X, Y, Z and the coordinates 
x, y, z by the same common factor, the entries of the matrix rescale with the same power 
and thus the eigenvalues scale). 

Let us solve the problem of the spectrum first in the asymptotic regime, where lets say 
(x, y,z) — > oo along a determined direction in M 3 keeping the X, Y, Z matrices fixed. By 
convenience, we can use rotation invariance to take z — > oo keeping x, y equal to zero. 

Then we have that 



We can compute the eigenvalues of H e ff by considering it as a perturbation theory of 
H e ff ~ — za z . The eigenvalues of this matrix are degenerate. There are N eigenvalues 
of values +z and N eigenvalues of value —z. These are very large. Since the spectrum 
is degenerate, to first order we need to resolve the splitting among the degenerate subset. 
This is done by looking at the perturbation terms in H e ff that commute with a z . The term 
that does that is Za z itself. So the transformation that diagonalizes H e jf along the two 
degenerate subsets is the same transformation that diagonalizes Z. 
We find that the leading order spectrum is given by 




(12) 



Etg(H eff ) = ±(z-\f) + 0(l/z) 



(13) 
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where Af are the eigenvalues of Z . The extra corrections of order 1/z are from perturbation 
theory: they result from 'energy denominators' and involve the components of X, Y. 

We find the familiar theme that the eigenvalues of the matrices X, Y, Z describe the 
positions of objects (distances) as viewed from infinity. Since the eigenvalues are continu- 
ous functions of the matrices, we find that the notion of distance by taking the minimum 
eigenvalue is a continuous function of the position. 

We are now interested in asking what happens when we are at distance zero from a 
configuration. 

This can happen in two ways: an eigenvalue of H e ff crosses zero, or the eigenvalues just 
graces zero and keeps its sign. A really interesting question is whether the spectrum of 
H e ff always has paired eigenvalues: if eigenvalues cross zero, this is not so. If eigenvalues 
are always paired, then every time one eigenvalue reaches zero from positive values, then 
another one reaches zero from negative values. We will define an index that counts possible 
crossings of zero from infinity. At infinity, the spectrum of H e ff is paired into positive and 
negative eigenvalues and to first order in perturbation theory they are equal to each other 
up to sign, obviously this implies that they both have the same number of eigenvalues. If 
an eigenvalue goes from positive to negative, the number of positive eigenvalues decreases 
by one, and the number of negative eigenvalues increases by one. Similarly in the other 
direction. We want the index to be zero at infinity and to change by one by each crossing. 
The definition of the index is given by 

I((x,y J z)) xx , z = r ^^ (14) 

where n + the number of positive eigenvalues of H e ff, and n_ is the dimension of the space of 
negative eigenvalues of H e ff. The index is a locally constant function (after all, eigenvalues 
of matrices are continuous functions of the entries) that can only change values at locations 
where H e ff has null eigenvalues. If for a configuration we have that I(x,y,z) ^ 0, we 
know that on any path connecting x, y, z to infinity there are crossings of zero and thus the 
location (x, y, z) is surrounded by the noncommutative object characterized by X, Y, Z. 

Such an index was defined in [21] for the position x,y,z = 0. It was called a Bott 
index. In their formulation, they were dealing with approximations to a sphere, where 
X 2 + Y 2 + Z 2 ~ 1 and the introduction of Pauli matrices was an auxiliary construction 
in mathematics. The matrices X, Y, Z represented observables in a quantum system where 
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only finitely many states are allowed and hence position observables become finite matrices. 
They were also restricted to have small commutators. The operator i? e // i n that case would 
square to something that was very closed to the identity so all eigenvalues of H e ff would 
need to be very close to ±1. Counting positive and negative eigenvalues is an invariant under 
small deformations that prevent the eigenvalues from getting too far from ±1. The index 
as interpreted in that case was an obstruction to localizing the states on a sphere (making 
X, Y, Z strictly commute), by demanding that \\X 2 + Y 2 + Z 2 — r 2 \\ < 5 by deformations 
of X, Y, Z that keep this property and a bound on their commutators is implemented. The 
spectrum of the operator X ■ a is also used in numerical studies of noncommutative field 
theories ( see [23] for a recent example), and one can also use the operator X ■ o to define a 
fuzzy sphere by studying a single matrix model of IX x 2N matrices with a constraint [24J . 

In the case we have described here the index is dictated by the dynamics of fermionic 
degrees of freedom on D-branes. There are also no restrictions on the size of commutators. 
These ideas can be extended further to higher dimensions and matrices with various restric- 
tions following the ideas in [22]. Such a generalization is beyond the scope of the present 
paper. 

Here are basic properties of the Index function (some of these already appear in the work 

my- 

1. The index is an integer. At infinity the index vanishes ( as we computed already). 
The index changes by ±1 if a single eigenvalue crosses zero. It changes by integers if 
many eigenvalues cross zero. 

2. Orientation: the index defines a collection of oriented closed surfaces. The surfaces 
are the locus where the index changes value. The orientation is defined by going from 
larger values to smaller values of the index (this includes the sign, thus —1 > —2 
etc). The surface set itself is obtained from the zero locus of a polynomial in (x,y,z) 
obtained by taking determinants. These surfaces will be called membranes or D-branes 
interchangeably. 

3. Additive property. Given two configurations Xi,Y\, Z\ and X2, Y2, Z2, we can consider 
a new configuration given by taking direct sums X 3 — X\ © X 2 , Y 3 — Y\ © Y 2 , Z 3 — 
Z\ © Z 2 . The index is additive under such constructions 

y, z))x 3 ,y 3 ,z 3 = I((x, y, z)) Xl , Yl , Zl + I((x, y, z))x 2 , Y 2,z 2 (15) 
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and the set of surfaces with orientation is also additive under this operation. 

4. Orientation reversal. This states that we can reverse the orientation of a surface 
without affecting its shape. This is done by considering the complex conjugate to the 
matrices X, Y, Z. In equations, we have that 

I((x, y, z)) x ,y,z = -I((x, y, z)) x *,y*,z* (16) 

This property is less obvious. A proof is as follows: A matrix and its transpose have 
the same eigenvalues. Thus (X — x) <g) a has the same eigenvalues as (X — x) T cg> a T . 
Now, X T = X*, so we can substitute. However, for Pauli matrices we have that 
a T ~ —a after a unitary transformation in the spin one half subspace. Thus, we have 

— * — * 

that the eigenvalues of (X — x) <g> a are equal to the eigenvalues of (X* — x) <g) (—<?). 

— * — * 

This is, the matrix (X* — x) <8> a has the same eigenvalues as [X — x) <8> o but with 
signs changed. This exchanges n + and n_ and reverses the index. 

5. If X, Y, Z are real, then I((x, y, z)) = everywhere. This is a corollary of the orienta- 
tion reversal property. Obviously for such configurations we have that I((x, y, z))x,y,z = 
—I((x,y,z))x*,Y*,z* = — I(( x i Vi z ))x,y,z- From which the result follows. This in par- 
ticular holds for collections of zero branes: direct sums of one dimensional problems. 

6. The index is covariant under rotations, translations and dilatations of the system. 
This follows from the similar properties that H e ff has. 

7. The index is not trivial: there are matrix configurations (X, Y, Z) for which I((x, y, z))x,y, 
0. We will explore these in the next section. 

If we instead work with the BMN matrix model we get an effective Hamiltonian given by 

3 

H eff = -za z + (Za z + Xa x + Ya y ) + -crW (17) 

The extra term causes trouble with scaling the surfaces, and with being able to change the 
sign of the eigenvalues by complex conjugation. This way various of the properties above 
are broken. For example the change of orientation does not happen automatically, and 
the corresponding Index does not behave as nicely. We still get translation and rotation 
covariance. The index still vanishes when a probe is at infinity, but one can check that even 
for 1 x 1 matrices, the index changes when the probe which is used to define the index is 



17 



on top of the 0-brane that described the configuration. This is the Myers effect in action 
|2Uj . Indeed, as far as fermions are concerned, the presence of a background RR flux changes 
the Dirac equation, and an example computed by one of the authors of the paper can be 
found in [25]. In that example the displacement of the location of the fermion zero modes 
was required in order for configurations to form tori that were BPS. In the present case, the 
structure of the gamma matrices follows the background flux in the BMN model [7]. 

On the other hand, in this case many fuzzy spheres are ground states of the system and 
one expects that these solutions survive as time independent configurations. Also, many of 
these can be made to oscillate slightly so the membranes can persist indefinitely. 

IV. FUZZY SPHERES AND EMERGENT SURFACES 
A. Fuzzy spheres 

Now that we have defined an index, let us consider some special examples of the index 
computation. We will start with a fuzzy sphere and ask about the index at the center of 
the sphere. The fuzzy sphere is defined as follows. Let L 1 ' 2,3 be the angular momentum 
matrices of the irreducible representation of SU{2) of spin j. These satisfy the identities 

[L\ U] = ie i]k L k (18) 

The maximum eigenvalue of L 3 are ±j. Consider the following set of 3 matrices built by 
the following combinations: 

rr> qr» qr» 

X = -L 1 , Y = -L 2 , Z = -L 3 (19) 
J J J 

This is called a fuzzy sphere. The maximum eigenvalue of Z is \r\. Thus one could argue 
that the sphere has radius |r| (as seen from infinity as in our large distance computation 
in the previous section). Notice that X 2 + Y 2 + Z 2 = ^±^r 2 . Thus one could also argue 
that the radius of the sphere is given by f = + j)|r|. These two become identical in 
the large j limit, but at finite j there is some discrepancy. However, it is natural to believe 
that there is a well defined surface near this radius that surrounds the origin and that is our 
candidate for a locus where an eigenvalue changes sign. 

Let us prove this assertion by computing the index in the center of the configuration, at 
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x = y = z = 0. The effective Hamiltonian we have to deal with is then given by 

Heff = ~X-B (20) 

This is the same type of problem that shows up in spin-orbit coupling in the hydrogen atom. 
The important thing is that this is spherically symmetric, so it makes sense to decompose 
the Hilbert space Hilbu g into irreducible representations of SU(2). The big Hilbert space is 
given by 

Hilb big ~ (j) ® (1) ~ (j + 1) © (j - 1/2) (21) 

and it decomposes into two irreducible representations of 577(2). For each of them, we have 
a common eigenvalue of H e ff. Moreover, i7 e // is traceless. This can be proved in general 
because the Pauli matrices themselves are traceless. Thus, the two possible eigenvalues of 
H e ff have the opposite sign. One is positive, and the other is negative. This depends on 
the sign of r. Let us choose the sign of r so that n + > n_. The number of eigenvalues of 
the bigger representation of SU(2) is n + = 2j + 2, while those of the smaller representation 
are n_ = 2j. These are the dimensions of the two irreducible representations of 577(2) 
appearing in the tensor product. We obtain that 

'((0, 0, 0)) Fuzzy Sphere = = 1 (22) 

We already knew that the index was an integer, and that the typical change should be by 
±1. Here we find an explicit example where the index changed by one somewhere between 
the origin and infinity. Because of spherical symmetry, the index changes value at a fixed 
sphere radius. A direct computation carried in the appendix shows that the radius at which 
it happens is given exactly by |r|. We thus find that the radius is governed by the maximum 
eigenvalue, rather than by the value of X 2 + Y 2 + Z 2 . Indeed, if we use the definition of 
distance from the origin that is obtained from the spectrum of i7 e // we find that the distance 
is equal to |r|. Indeed, with the spectral definition of distance we used, we find that the 
distance from any point in space to the sphere is the one that is obtained by elementary 
geometry. 

Obviously, we can also set up direct sums of concentric fuzzy sphere configurations of 
various radii, so we can get configurations where the index is arbitrarily large. For such 
configurations the index counts the (minimal) number of sphere layers that need to be 
crossed to get out of the center. Since the index counts with sign, surfaces (which we call 
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membranes) of different orientations can be present and the index itself represents a lower 
bound on the number of layers that need to be crossed. 



B. Prom a sphere to a torus 

Here we detail how to make configurations that lead to a fuzzy torus embedded in three 
dimensions. The idea is to begin with a fuzzy sphere and to deform the matrices in a simple 
form to go from a sphere to a torus. The basic idea is to follow the construction of the giant 
torus as described in [2S] (other examples of embeddings of Riemann surfaces in IR 3 can be 
found in [27], and in [2S] one can also find a different example that interpolates between 
sphere and tori). In the case of the giant torus, one is supposed to add strings with maximal 
angular momentum to a sphere until the geometry transitions to a torus. To do this, it is 
convenient to use matrices defined by 



X ± = X ± iY (23) 



and in the other direction 



A' = ^ (24) 

The matrices X ± in the fuzzy sphere of spin j case are rescaled ladder operators for spherical 
harmonics. X are adjoints of each other. In a natural basis for a sphere, we have that 

X+ = r ^(j(j + i)- a (a + l)6 bia+1 (26) 

The labels a, b go from ]■■■—]■ 

In the matrix X + , if we quantize fluctuations of the fuzzy sphere (see for example [29J), 
one can check that the different diagonals of the matrix carry different amounts of angular 
momentum in the z direction. They differ by one unit, and the diagonal where X + has 
entries carries no angular momentum in the z direction. When we condense various of these 
fluctuations, we simply replace them by an expectation value which becomes just a number 
multiplying the appropriate fuzzy tensor harmonic. Since we are looking to maximize the 
angular momentum of the fluctuations, the deformation we seek is given by 

X+ = r ^(j(j + l)-a(a + l)5 b , a+1 + r^5 bJ 5 a ^ (27) 
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and we then take X~ = (X + y. We are using the index convention for the matrices that is 
associated to the L z spin of the SU(2) representation of spherical harmonics. The self-adjoint 
matrices X,Y are built from the same linear combinations as above, after the deformation. 
The parameter r just rescales the full solution, so we can ignore it. The parameter (3 then 
controls the geometry. For (3 = we have a sphere. Indeed, the topology of the sphere is 
preserved for some values of around zero. We have seen numerically that the topology 
changes at the precise value (3 = j, this is not essential for our discussion. 

Another thing to notice is that the presence of (3 breaks the rotational symmetry to Z 2 j+i 
which is the rank of the matrices. This can be understood from the spin of the excitations 
around the z axis: it is the unbroken symmetry associated to condensing the state with 
highest spin along the z axis. Thus the torus shape is not invariant under full rotations 
along the X, Y plane. The simplest case where the family of surfaces we get seems to 
contain a torus is for 4x4 matrices. A figure for the case of 6 x 6 matrices is presented in 
figure [Tj 




FIG. 1. A fuzzy torus, for r = j = 5/2, (3 = 2.55. The symmetry is easily visible. 

We should also notice that in our case it is obvious we have a torus. In other setups, to 
argue for the genus of the surface is more involved, using an approximation to Morse theory 
on the surface [30], and the result is inherently more fuzzy, or explicitly requires taking a 
limit of large matrices |32j. One can also obtain more standard fuzzy tori as zero energy 
configurations in higher dimensions by studying beta deformed matrix models [31J. 
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C. D2-branes 

The main characteristic of D-branes is that they carry a connection on their worldvolume. 
This is the familiar statement that the open string sector has a massless spin one particle 
on the D-brane worldvolume. Here, we want to show that the geometry that we deduced 
for these surfaces carries the information of a line bundle on it. 

The idea is rather simple. The surfaces are defined by the vanishing of an eigenvalue 
of H e ff, as calculated for equation (10). Obviously, for a single zero eigenvalue there is a 



corresponding eigenvector. Let us call it ip (x,y, z). 

The eigenvector if> , normalized to unity, is well defined up to a U(l) phase. This is the 
familiar symmetry for states in a Hilbert space in quantum mechanics: a global phase for 
the full wave function is not measurable, as physical states are rays in the Hilbert space. 

Now, for each position on the surface there is such an eigenvector. This changes con- 
tinuously when we vary the position along the surface (x,y,z), as the eigenvectors are also 
smooth functions of the matrix entries. One can easily understand this fact by the fact that 
the eigenvectors can be calculated using perturbation theory in quantum mechanics. 

One can construct a bundle from these ipo(x, V, z). One defines sections of the bundle by 
functions multiplying ^0(^5 V, z )- Because the phase of ip is ambiguous, we have to choose a 
phase by patches on the surface, and between patches there are transformation rules for ip Q . 

One can also define a connection on the patches. This is done by the familiar Berry 
phase, defined by 

= -iv^*{x,y,z)d fl M^y,^) = -iv^iAWo) (28) 

where is a tangent vector to the surface. Obviously, this defines the connection of a line 
bundle on the worldsheet. Thus, at least in principle, the membrane behaves exactly like 
we would expect a D2-brane to behave. At this stage, it is not clear the Berry connection 
that one would compute this way is just the connection that open strings feel, or if this is 
further twisted by the tangent bundle on the surfaces that were defined as we prescribed. 

The full exploration of the curvature on these bundles and the precise connection to D- 
branes is beyond the scope of the present paper. We will show later that there is further 
evidence for physical states feeling a connection on the membrane worldsheet when we 
intersect two of these objects. 
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V. A LINKING NUMBER 

Now that we have defined a geometric object for a collection of 3 hermitian matrices, we 
can do something more. We can take two such objects and ask how they are related to each 
other. Indeed, in the matrix model setups, each of them would be a configuration of branes, 
so the spectrum of strings stretching between them becomes interesting from a dynamical 
point of view. One can define a linking number that for zero branes at a position x reduces 
to the index we defined in previous sections. 

The idea is to take the matrices X 1: Y 1: Z x of rank r\ and X 2 , Y 2 , Z 2 of rank r 2 and define a 
matrix analog the Hamiltonian H e ff(x,y, z)x 1 ,y l ,z 1 , where we replace (x,y,z) by hermitian 
matrices (X 2 , Y 2 , Z 2 ). If the matrices commute with one another, we want the H e ff operator 
to give us an operator that acts as if e // on the direct sum over the eigenvalues of X 2 ,Y 2 , Z 2 . 
One easily sees that the following Hamiltonian does that: 

H^f f (X 1: X 2 ) = {X 1 <g> l r2 - l n ® X 2 ) <g> a x + (y x) + (z x) (29) 

Notice that once is defined this way, it does not matter anymore that the X 2 , Y 2 , Z 2 
matrices commute with each other. 

Then, the definition of our linking number is given by 

LW[(X 1 ,Y 1 ,Z 1 ),(X 2 ,Y 2 ,Z 2 )] = (30) 

It's easy to prove that is antisymmetric in the entries. This is because tensor product 
spaces A <g) B are equivalent to B <g> A as Hilbert spaces. If we think of these spaces in 
tensor notation, the equivalence is a reordering of the indices. The Hamiltonian then 
changes sign (more precisely, H^(X , X') is unitarily equivalent to —H^(X',X)) when we 
exchange the triples X 1 and X 2 . 

There is a second linking number that one can define, by changing a brane by an antibrane, 
this is, reversing orientation: 

H^f f = (Xi ® l r2 - l n ® X*) ®a x + {y^x) + (z^x) (31) 

(2) 

Again, if X 2 ,Y 2 , Z 2 commute with each other and are diagonal, we can not distinguish H e ^ 
from H^j. But if the matrices do not commute with each other, we can. The definition of 
the second linking number is 

L^[(X 1 , Y 1 , Zi), (X 2 , Y 2 , Z 2 )\ = (32) 
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This is symmetric in the exchange of (Xi,Y]_, Z\) and (X2,Y2, Z-i). This uses the antisym- 
metry of combined with the change in sign of the index upon complex conjugation 
discussed in previous sections. It turns out that when considering the dynamics of fermions 

(2) 

as given in the BFSS matrix model, it is the spectrum of H\tt that controls the physics 
|33j . This is because the matrix multiplication rules on commutators translate to needing 
to take the transpose of the matrices X 2 ,Y 2 , ^2, which is equivalent to using their complex 
conjugates. This is also equivalent to saying that the fermions transform as a fundamental 
under one set of branes and an antifundamental with respect to the other set of branes. 

Also notice that if we move one of the objects and take them to infinity (by adding 
multiples of the identity matrix), then at infinity both of the linking numbers are zero. 
Also, if we shrink one object until it is point-like (by making X 2 , Y 2 , Z 2 proportional to the 
identity matrix, with coefficients x 2 , 1/2, z 2 ), then the linking number is r 2 times the index 

I(x2,y2,Z 2 )x 1 ,Y 1 ,Z 1 - 

Also, one can use these same Hamiltonians i^// an d -^e// t° define a spectral distance 
between two such configurations, again by taking the eigenvalues closest to zero and taking 
absolute values. For infinite membranes touching each other in the IKKT matrix model one 
finds zero modes |34j . The effective Hamiltonian for fermions in that case takes a similar 
form to the BFSS matrix model. This is just as expected from the mode spectrum of brane 
intersections at angles [35]. When the intersections are extended and compact, the low lying 
modes at the intersection need to be quantized carefully and zero modes are not guaranteed. 
One would expect that the spectral distance then gives an upper bound for a geometric 
distance between the brane configurations. 

We will now give an application of the linking numbers. We will show that the linking 
numbers actually take into account that the surfaces that are defined by previous sections 
actually carry a connection for a line bundle on them that couples to physical states. This 
provides further evidence that the surfaces are actually behaving as D2-branes. This is 
easiest to check from the calculations in the appendix. 

The idea is as follows: take two fuzzy spheres and displace them relative to each other. 
For simplicity, we have them normalized so that the radius is equal to j and f , the spin of 
the corresponding representations of SU(2). This is natural in the BMN model for ground 
states. Let the displacement between the fuzzy sphere centers be characterized by b. Because 
of the high amount of symmetry, one can actually compute the index analytically and follow 
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the crossings of zero of the fermion eigenvalues in a lot of detail. If the displacement is 
b, along the z axis, and the fermions are decomposed into fuzzy spherical harmonics with 
respect to both fuzzy spheres, one finds that the eigenvalues cross zero sequentially when 
b = j — f + £, where I is an integer between and 2j' inclusive. This means the index 
of the configuration where the big fuzzy sphere surrounds the smallest is exactly equal to 
2f + 1, which is the dimension of the representation of the spin j' set of matrices. This is 
expected: the small fuzzy sphere is made of 2j' + 1 DO branes, so that when they are all 
inside the big fuzzy sphere, we expect the index to be 2f + 1 times the index of the smallest 
representation. 

The first zero mode appears when the spheres touch each other for the first time, at 
displacement b = j — j'. As b advances further, the two fuzzy spheres touch each other 
along a circle. We expect the lightest fermions to be localized in this circle. So the problem 
effectively reduces to a one dimensional problem. As can be seen from the results of the 
appendix, the fermion modes with maximal angular momentum in each SU(2) representation 
of fuzzy spherical harmonics do not mix with other states, and their frequencies are given 
exactly by 

j-f + l-b (33) 

where I is an integer. These states are evenly split in energy, creating effectively a Kaluza- 
Klein tower of finitely many states (this is very similar to the Kaluza-Klein tower of tachyons 
between such spheres computed in the BMN model for such crossings in [E]). Such a Kaluza- 
Klein tower is an approximation to a quantum field theory for zero mass fermions on a circle 
(either leftmoving or right-moving depending on if the energy of the mode is positive or 
negative) in the presence of an holonomy around the circle (this can be translated to quasi- 
periodic boundary conditions on the fermions if we want to). The fermions can have zero 
eigenvalues if the holonomy is a multiple of 2n. This can be removed by a large gauge 
transformation redefining the notion of momentum on said circle. The important thing to 
notice is that if the corresponding surfaces that are intersecting have the properties of D- 
branes, in that each carries a connection A\, A2 , then the fermions that stretch between 
them feel the connection A\ — A%. Because of spherical symmetry, this connection can be 
computed from Gauss law by calculating the area of the sphere that the circle where the 
fermions lie enclose § Aids = f s Fida. The area of a sphere slice is proportional to the 
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height of the slice, hence f s F\da oc A which is linear in the vertical height of the slice. 




FIG. 2. Illustration of two intersecting spheres. The net connection seen by the fermions stretching 
between can be computed by calculating the net flux through the solid lines. 

Hence, shifts of 2n in the holonomy are equally spaced in b. Indeed, if the sphere is 
made of n DO branes, we expect the total flux through the sphere of this bundle to be equal 
to n = 2j + 1. However, we can also expect a curvature correction. If we think of the 
matrices as describing a lowest Landau level of endpoints on each sphere, in order to have 
n states we need a monopole flux equal to n — 1 (this is if the endpoints are treated as 
monopole spherical harmonics). This extra one is the contribution of the curvature of the 
sphere. One can check this way that the net flux for this connection through each sphere is 
2j and 2f respectively, as opposed to 2j + 1 and 2f + 1. Thus the net flux through a slice 
is proportional to the area, which is 2irj(2j — t) where t is the height of the slice. Since the 
total flux through each sphere is 2j, and the area is Airj 2 for each sphere, we get that the flux 
per unit height is constant and the same for both spheres. Thus, the flux for each sphere 
is linear in height with the same coefficient. This can be visualized in figure [2} As seen in 
the figure, the net flux that we need to compute is the one associated to the surface that 
has not been dashed in the graphic. This is proportional to j + j' + b as a function of the 
displacement. We need this number to be an integer in order to get the correct holonomy. 
We see then that the geometric argument matches the matrix computation. Obviously this 
is a simplified computation where it so happens that the flux per unit height on each sphere 
is the same. 
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Thus, the setup shows that indeed the surfaces we make are compatible with the idea of 
having a curvature of a line bundle on them for physical states that thread between them. 
There is another way to think about this that we already discussed: on the locus of positions 
where an eigenvalue vanishes there is a preferred fermion wavefunction for the hamiltonian 
H e ff. this is the zero eigenvector itself. This is only well defined up to a phase. If we want 
to patch these together to form a vector bundle over the surface, we need a line bundle 
connection so that this phase ambiguity is resolved on parallel transport. This is the generic 
case, but we can set it up so that the null eigenspaces are degenerate (thereby giving us 
multiple branes on top of each other). Thus, in general one will need a bundle connection 
to resolve these issues. Since the structure that we are analyzing involves the symmetries of 
a Hilbert space under change of basis, the connection in general will be U(n) valued for n 
coinciding branes. 

The last thing that we will do in this section is to give a more physical interpretation 
of these zero modes. The main idea, which we have hinted at already when we defined the 
index function, is that a crossing by zero represents a raising operator becoming a lowering 
operator and viceversa (for the particle conjugate). If we follow a ground state continuously 
past this change, the ground state is defined by a\0) = = b\0), where a is the lowering 
operator for particles (the ones with positive frequency), and b is the lowering operator for 
the antiparticles (the ones associated to negative frequencies in H e jf). After the crossing 
by zero, the state that follows by continuity of |0) is not ground state anymore. Instead, 
one of the lowering operators, let us say a a , becomes a (a raising operator). The state on 
the other side of the barrier will have a non-zero occupation number for a single fermion in 
the Hilbert space. This is, on crossing a zero, a fermion is created. This is nothing but the 
Hanany-Witten effect (and its various generalizations discussed in [SUH])- 

The linking number we defined then encodes the number of strings that are created (with 
orientation) when separating two objects that are partially inside each other. Or the number 
of strings that were created on bringing the objects together from infinity when they cross 
each other. This is done by following a vacuum adiabatically until exactly the point where 
the transition happens (where there is a degeneracy of vacua), and then following the state 
that is created after a crossing and which is not a vacuum any longer adiabatically as well, 
until further crossings where fermion zero modes determine a degeneracy of a Fock space of 
fermions at each level. 
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Notice that this interpretation in terms of the Hanany-Witten effect explains why the 
geometry is so sharp. The presence or not of strings connecting the two surfaces is easy 
to test: we check if we the fermionic ground state is gauge invariant or not. The Hanany- 
Witten effect has the property that the fermionic ground state is not always gauge invariant, 
so the presence of the strings is protected by topology. 

VI. ASPECTS OF MATRIX BLACK HOLES 

Recently, various simulations have been carried in matrix models to understand various 
aspects of the dynamics of black holes in holographic setups. The main idea so far has 
been to compare the numerical simulations in the BFSS matrix model with black holes as 
described in [36]. The numerical approach was initiated in the works [37J EH] and a lot of 
the thermodynamic static properties of the black holes have been matched in the quantum 
mechanics. The most impressive such agreement is in [32] • The BFSS matrix model has an 
infinite moduli space of vacua, so the thermal ensemble of these models is not well defined. 
This gives such calculations a systematic error. To have a better setup one wants a matrix 
theory with a well defined ensemble, and the BMN matrix model fits the bill. Numerical 
simulations using lattice techniques were carried out in |40j. Also, classical simulations of 
the dynamical evolution of the BMN matrix model have been carried out in [5] . 

All of these calculations in general give us a huge number of sets of matrices about which 
we can now start asking very geometric questions, for example: how many membranes are 
inside the black hole? One good reason to do this is that general consideration of black 
hole entropy for non-extremal black holes suggests that they are made of a gas of brane 
anti-brane pairs [12] and their excitations. 

Since our construction permits us to study the geometry of the typical matrices in a ther- 
mal ensemble, we can ask how does the black hole look like in the matrix variables. We will 
show an example of this based on data obtained from simulations similar to those reported 
in [5] , where we truncate to only three matrices as described previously. The simplest thing 
to do in order to understand the data is to compute the spectrum of H e ff(x, y, z)x,y,z given 
a collection of three matrices X, Y, Z that are obtained from such simulations. We show the 
results by fixing the matrix, and setting y, z = 0. A typical such result is shown in figure [3} 

What we should notice is that there are various crossings of zero, and that there is a 
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Eigenvalues H^j- 




FIG. 3. Eigenvalues of H e ff(x, 0, 0)x,y,z from a typical configuration of matrices after thermaliza- 
tion, when varying x. The matrices have rank 21. 

region in the center of the configurations where the eigenvalues of H e ff do not seem to have 
a gap in them. That is the black hole region of the configuration. For large x we see that the 
eigenvalues behave as parallel lines and this matches our expectations based on perturbation 



theory from equation (13). We clearly see various crossings of zero, mostly because of the 



shift of flux. Also, in the region without a gap the eigenvalue distribution appears to have 
a well defined density of eigenvalues. 

One can show that if one varies the classical temperature of the BMN ensemble T and 
one makes it large (so that the quartic term in the potential dominates over the cubic and 
quadratic terms), then the matrices themselves grow roug hlylikeT 1 / 4 . This is an application 
of the virial theorem applied to the BMN Hamiltonian. Similarly, one can show that at T 
fixed, and varying N, the matrices grow only like iV 1 / 4 . The virial theorem would just state 
that 

E kin ~ N 2 T ~ E pot ~ Tr([X, X] 2 ) ~ Tr(X 4 ) ~ Nx A (34) 

this is, the kinetic energy is proportional to the potential energy, which is roughly character- 
ized by the typical eigenvalues of the matrix X, which we call x in the equation. Since there 
are N such eigenvalues and the potential is quartic, we expect that the result is roughly 
Nx\ We then get x ~ N^T 1 / 4 . 

This means that the matrix H e ff also grows roughly like iV 1//4 , so it's eigenvalues scale like 
iV 1//4 . If we assume that in the ungapped region the Hamiltonian H e ff behaves like a random 
matrix in that it has a well defined density of eigenvalues when we take N — > oo, then since 
we have 2N eigenvalues, the eigenvalue density near the black hole region near zero grows 
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like iV 3//4 . Because of the flux contribution to the fermion mass (the Myers effect [20]), the 
center of the configuration is displaced from zero: the eigenvalue configuration is centered at 
3/4 in our units (see appendix). Far away, half the eigenvalues are above the x axis, and half 
are below, so roughly p(0) ~ 0(iV 3//4 ) eigenvalues cross zero. This means the inside of the 
black hole is full of branes that have been polarized, all with the same orientation. As the 
temperature is increased, the eigenvalue distribution becomes wider and fewer eigenvalues 
cross zero. The N dependence is still correct, but there is also a temperature dependence. 
This polarization into D-branes is mostly because of the Myers effect. Remember that in 
this problem we have truncated to three matrices, so we are actually working on an orbifold: 
we are only using H e ff with only Pauli matrices. The true fermion Hamiltonian in the BMN 
matrix model uses all matrices and will have different characteristics. Thus, if we truncate 
this way we are working with something that resembles more a brane-world black hole (we 
can't move it away from some locus). 

A. Black holes absorb matter 

An important characteristic of black holes is that if one throws matter at them, then the 
matter does not come out at the other side. Let us throw a fractional DO brane at such a 
(orbifold) black hole. Notice that we need to do so in the orbifold of the BMN geometry. 
We can require the fractional DO brane to be at a large distance from the black hole (let us 
say k times the size of the black hole itself). The energy of such a DO brane in the BMN 
matrix model that starts at rest is of order k 2 ~N x l 2 . To estimate this we just look at the 
quadratic potential term for the X 1 ' 2,3 matrices. 

If we throw a fractional DO brane to the black hole as described above, at each zero 
eigenvalue crossing a string is created due to the Hanany- Witten effect we have discussed 
previously. This is identical to the creation of strings observed in the one dimensional model 
for D-brane scattering studied in [JT] • There are about N 3 ^ 4, such strings per DO brane that 
are created (as many going in as out). When we reach the end of the matrix configuration, 
these strings have a length of order N 1 ^ 4 , so the energy stored in these strings is of order 
hN. So long as HN >> k 2 N l l 2 , we find that the DO brane does not have enough energy 
to come out: it gets transfered to the strings. Notice that this depends on h. When we 
take N — > oo, it's clear that the strings win over the initial energy of the DO brane. Thus 
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we find that the matrix thermal configuration does become a very good absorber: every 
fractional DO brane that is thrown at it is eaten. Indeed, if we throw bigger objects at the 
thermal ensemble (let us say made of M fractional DO branes), the effect is proportional to 
the number of fractional DO branes making the object: the linking number will guarantee 
that. Thus all objects are absorbed with the same efficiency This is very similar to how 
black holes actually operate. 

Notice also that usually in these dynamical setups if fermions are created by dynamics, 
then so are bosons. The accounting might be different, but they usually follow each other 
somewhat. Indeed, in the BMN model alone, the presence of tachyons in some regions of the 
dynamics can generate large numbers of bosons [33]. Thus one should also expect bosonic 
modes to be created by dynamical mechanisms (the modes become non-adiabatic) rather 
than by a simple topological argument in general. 

This simple accounting of how objects are absorbed that we found is different than other 
approaches that presume the formation of a tachyon in an ensemble [33]. Maybe an effective 
tachyon can be thought of as a collective effect of all these fermions and bosons. 

We should remind the reader that we should not take the arguments above based on 
generalizations of the Hanany-Witten effect very seriously for the full BMN matrix model. 
There the dynamics of the other matrices might change the physics substantially, as there 
we expect these D2-branes to fluctuate in the transverse directions. Thus, the topology of 
the Hanany-Witten effect would only be available for D8 branes, rather than D2-branes, so 
that there is no background flux reason to polarize D8 branes in large numbers. 

VII. CONCLUSION 

The 16 supersymmetries of the BFSS matrix model can be reduced down to 3 supersym- 
metries, removing six of the nine bosonic matrices and thereby giving three matrices which 
capture some dynamics of the full theory. Such a reduction can be obtained by orbifolding 
six of the nine directions. If the orbifold is chosen with respect to a action, then chiral 
fermions arise. Chiral fermions give rise to anomalies in four dimensions, and it follows that 
the fermions can encode some topological information in the reduced matrix model. This 
information can be used to study the geometry of membranes formed by thermalized black 
holes in numerical simulations of the matrix models. 
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Appending a DO-brane probe, described by a point in M 3 , to the three relevant matrices of 
the matrix models allows us to look at the dynamics of the theory. An effective Hamiltonian 
derived from the interaction between chiral fermions, described as fractional branes different 
than those we're probing, and the branes in the configuration we are probing describes how 
fermionic strings are created between these branes and the DO-brane. The energy eigenvalues 
are proportional to the string length and thus the minimum eigenvalue can be interpreted 
as the minimum distance between the DO-brane and the noncommutative configuration. If 
an eigenvalue equals zero, then the probe is intersecting the membrane. Passing through the 
membrane changes the number of positive and negative eigenvalues by integer increments. 
An index function is built that captures these crossings by taking the difference of the count 
of positive versus negative eigenvalues. 

The index inherits the symmetries of the Hamiltonian from which it is derived; it is co- 
variant under rotations, translations, and dilatations and is gauge invariant. The continuity 
of the eigenvalue functions of a matrix imply that the index function is locally constant and 
defines closed oriented surfaces. It has been shown that the index is zero at infinity. The 
index is also additive amongst direct sums of different matrix configurations. Finally the 
index has a transformation that reverses the orientation of the membranes. These properties 
are all true in the BFSS matrix model. In the BMN matrix model, there is a mass con- 
tribution to the fermions caused by the dimensional reduction of the 9 dimensional gamma 
matrices. This ruins the orientation reversal and scaling properties of the index function. 
Furthermore, the added mass changes the shapes of the membranes slightly, which is directly 
related to the Myers effect [20J. 




Physically we can speak of the eigenvalues of the effective Hamiltonian as representing 
the energies of the fermions in different regions between the membranes; positive eigenval- 
ues corresponding to fermion creation and negative eigenvalues corresponding to fermion 
annihilation. As one crosses a membrane from higher to lower index, a fermion creation 
operator is transformed into a fermion annihilation operator. In the BFSS matrix model we 
can view this as a generalization of the Hanany-Witten effect. In the BMN matrix model the 
dynamics of flux changes the results and many of the properties of the index are modified. 
This can be ascribed to the Myers effect [20]. These crossings define the surfaces in M. 3 . We 
showed configurations that correspond to both spheres and tori. We were also able to show 
that the surfaces carry the information of a line bundle on them with connection (which 
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can be calculated using Berry phase arguments). This shows that the membranes we found 
really behave like D2-branes. Our exploration of this issue was very sketchy, so finding how 
to make this correspondence precise requires more work. Indeed, we would need to see if 
the connection we computed also includes information of the tangent bundle of the surface 
or not and how to separate that part from the D-brane worldvolume spin one excitations. 

We were also able to generalize this index to a linking number between two such config- 
urations. The linking number is also interpreted in terms of the Hanany-Witen effect: it 
counts how many such strings are created when trying to separate the two configurations 
away from each other. 

Finally we were able to show that these surfaces can be used to analyze numerical data 
from simulations in the BMN matrix model and in more general setups and the data show 
that one can make contact with conjectures about the structure of black hole interiors as 
made from brane-antibrane systems. We were also able to show that with the Hanany-Witten 
effect, the fermions created on these surfaces could be used to stop a probe DO-brane particle 
in a simple model. Thus it is clear that these modes can give us a handle on black hole 
dynamics that do not require much effort. 

After solving the problem of embeddings into three dimensions, it would be interesting 
to understand how this same story plays out in higher dimensions. One of the ways in 
which extended objects are understood is in terms of Berry phase dynamics |45j ( for a 
more recent discussion see [46 J and references therein). The Berry phase can lead to a non- 
trivial vector bundle structure of the fermion excitations connecting a probe to a brane. One 
can expect that if topology requires that this structure becomes degenerate at various loci, 
that these loci describes extended objects: again, one has to look for fermion zero modes 
depending on position and at least in principle it should be possible to predict that there are 
degenerations in some setups. However, the story might be much more complicated, as we 
might require to have more than one fermion zero mode simultaneously to describe this locus. 
The Berry phase dynamics associated to that setup would then be non-abelian. It would 
be nice to understand this better. This might also require using extended probe branes 
to see the effects. The general question will then be to understand generalized versions 
of the Hamiltonian ([I]) and the general structure of degenerations. Our construction also 
suggests that in these general setups there can be a similar linking number so long as one 
can guarantee crossings of zero of the eigenvalues of H e ff. One can show that for even 
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dimensions (an even number of matrices) the spectrum of fermions starting from a DO brane 
probe to a configuration is mirrored: for every positive eigenvalue there is a negative one. 
This is because one can find a matrix similar to 7 5 in four dimensions that anti-commutes 
with H e ff as given by the generalization of equation 0. This suggests that ideally we 
should work with an odd number of matrices to make the existence of zero modes plausible 
for somewhat general configurations. 

A second thing that is interesting to study is how to recover the matrices given the 
surfaces (perhaps with additional information on them) and the rank of the matrices. One 
could also ask if the surfaces we obtained move in a way that closely resembles the membrane 
dynamics once we turn on the dynamics. This might be important to understand 1/N effects 
in matrix theory. Also, we found that in general we could reverse orientations of branes by 
using complex conjugation. It would be nice to understand if a brane-antibrane pair in 
these models generally leads to tachyons on their worldvolume and it would be interesting 
to analyze how tachyon condensation would progress in these setups. Also, it would be 
interesting to understand this issue with a probe DO-brane on top of a D2-brane: do we 
always get tachyons in this way? 

Also, the ideas found in [22] suggest various generalizations to different types of matrices. 
These ideas have applications in condensed matter physics and the connections we found 
with string theory ideas might provide interesting ways of analyzing the condensed matter 
systems and their dynamics. Such changes of the structure of matrices are natural when 
considering orientifolds. Thus our arguments should generalize to those setups. 
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Appendix A: BFSS and BMN model conventions 



The following was taken from I2H1 H3] using a mix of conventions: 



S — Sq + S m 
S = I dtTi 



(Al) 



^ 2(2/2) v uv ; 2 i 



j,fc=i 



9 1 

+ ^-(2i?)(*V[^>]) 



dtTr 
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J=4 



E f ^ 



(A3) 



The fermion representation we choose to work in is not explicitly real, and so \I/ T is replaced 
by \[d (see section ph. Rescale to get rid of R: 
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and again to get rid of /i: 



(A5) 
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S — So + S mass (A6) 

r 9 • 9 9 

f ^ ~ — 1 ~ J 1 « — 1 



i=i 



S, 



mass o 

g 2 



dtTr 



\j=i i=4 / 



+ 



i /3 



2 V4 



i,fe,i=i 



(A7) 
(A8) 



Notice that if we start with // = 0, then in the rescaling we just modify g and we find that 
the BFSS lagrangian has no free parameters (this is the statement that the gauge coupling 
in + 1 dimensions is dimensionful, so that there is no dimensionless coupling constant). 

In the A t = gauge, the covariant time derivatives become ordinary time derivatives. 
Relabel the (j>> by X 1 . Define X i = <p l for i = 1, 2, 3 and Y a = <p a for a = 1, . . . , 6. The 
bosonic action takes the form 



S 



B 



dtTr 



(X 1 ) 2 + (Y a f - (X 1 ) 2 - ^(Y a ) 2 - 2ie ijk X i XiX k - ^[X 7 ,X J ] 2 



(A9) 



The fermionic action becomes 



S F 



dtTr 



2 2 W 



(AlO) 



This is how the action is written in 



Appendix B: Fermion Decomposition 

This section comes from appendix A of reference Decompose the 16 component 

spinor as 

50(16) -> 50(6) g> 50(3) ~ 517(4) <g> 5f/(2) 
16 ^ (4® 2) © (4(8)2) 

tf-^Vla,^ (Bl) 

where 7, J are fundamental 5£/ (4) indices and a, (3 are fundamental 577(2) indices. The 
spinors obey the reality condition 
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which allow us to write the spinors in the stacked form 

* - ( *" ) (B3) 

The matrices g l3 are introduced to relate the inner product of SU(4) to the vector of SO (6) 
which satisfy 

g a {g ] ) b + g\g ] ) a = 25 ab (B4) 

The gamma matrices are then written as 

I -a* <g> 1 \ / 1 ® g a \ 

7* = , 7 a = . (B5) 

y a l ®I) \l®(g a y y 

The terms in the Lagrangian then decompose as 

!*t£> tf i^ Ia D i> Ia (B6) 

1*^123* -> ^N-ia (B7) 

i*V[^, *] -V> tJa *2f [**, ^] (B8) 

V 7 a [X a , vl/] -> Ie a ^ t/a ^j[y a , *P W ] - ^e^Vja^ 1 ) 01 - 7 ^ , W (B9) 
The fermionic part of the action (in the Aq = gauge) may then be written as 



Sf = I <It I 

2 



(BIO) 

Notice that the coupling to the X variables uses just the Pauli matrices after this de- 
composition. Also, a t/> spinor is always paired with its conjugate. If we perform orbifolds 
that are chiral, this structure remains, but the other mass terms that do not preserve four 
dimensional chirality might be eliminated. 

Appendix C: Fermionic Modes between displaced fuzzy spheres 

For our paper we want to consider computing fermionic modes between two fuzzy spheres 
in the BMN matrix model that have been displaced as described in [13]. We want to restrict 
to a chiral projection of the modes between two such fuzzy spheres. First we will setup some 
conventions for the fermionic modes of a single fuzzy sphere. Then we work with the more 
general problem. 
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1. Diagonal Fermionic modes 

This section is essentially a repeat of section 5.2 of [2H] using the conventions of 
The SU(4) indices are dropped as they do not come into play at all during the following 

calculation. We take the following conventions for the spherical harmonics and the angular 

momentum generators: 



Af = y/(l + m )(l-m + l) 
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where L = L\ ± %Li- We expand the fermions as 

lm 

The potential in the presence of the bosonic VEV's becomes 
V F - jTr(^a) = Tr [tf°o# [L\^]\ 
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= Tr 
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The eigenvalues of the matrix plus 3/4 give the mass spectrum. Note that —I < m,m' < I. 
The 5's tell us that this matrix has 21 two by two blocks and two one by one blocks where 
m,m' = I and m, m! = —I. Each of the one by one blocks yield the eigenvalue /. The 21 two 
by two blocks can be parametrized according to m from — I to I — 1. They are given by 

(C2) 
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The eigenvalues are given by the characteristic equation 

= (m - A)(-m - 1 - A) - (I - m)(l + m + 1) = (A - l)(X + I + 1) 

Thus the eigenvalues are / and — (I + 1). This means that the mass spectrum is M = 3/4 + 1 
with degeneracy of 21 + 2 and M — — (I + 1/4) with degeneracy 21. Note that there are two 
more positive eigenvalues than negative eigenvalues. 



2. Off-diagonal modes 



Here we follow the procedure of the previous section and that in |43]. Expand the off 
diagonal modes in fuzzy monopole harmonics 







to v (* m )n f 
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The following commutators are necessary 
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Substituting these expressions into the potential, taking the chiral projection and finally 
taking the trace we have 



m-b A l ™5 m i m A I 
A^m+i -(m-b) ) U 



)l 






, ^Im J 



(C4) 



Notice that this essentially produces the exact same matrix system as for the diagonal 
modes except with different diagonal elements. Also, half spin objects are allowed as when 
we decompose into the tensor product we can get different spins. There is a one by one 
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block with X = I — b, another with A = I + b, and 21 two by two blocks. The matrix for these 
blocks is 

(m-b A lrn \ 

(C5) 

l A i m _( m+ i 
with — I < m < I — 1. The eigenvalues satisfy 

= -(m - b - A)(m + 1 - b + A) - (/ - m)(l + m - 1) 
= A 2 + A - 1(1 + 1) - b{b - 1) + 2mb 

Solving for A gives 

A = -i ± v 7 ^ - m)(Z + m + 1) + (b - m - 1/2) 2 (C6) 
Thus the full modes are with — I < m < I — 1 (each with degeneracy two): 

m = i ± a/(Z - + m + l) + (6-m - 1/2) 2 

We also have two other modes corresponding to the one by one blocks of the mass matrix: 
m = 3/4 + /±6. These yield zero modes for the right value of b. That is we have zero modes 
when 

b = ±(1 + 3/4) 

for the modes with the greatest angular momentum in the z direction for a given value of 
£. These zero modes are correlated with the modes that become tachyonic for bosons in the 
same type of configurations found in [13]: they are objects of maximal spin fixing £. 
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